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cn '. Abstract 



are used to characterize the quark or gluon propagator, or other Green functions of higher 



K»" ! In the presence of the nontrivial QCD ground state or vacuum, nonlocal condensates 

m 

O . order. In this paper, we wish to show that, by taking the large Nc limit (with Nc the 

l> 

^ ! number of color) in treating higher-order condensates, a closed set of coupled differential 



Ph! equations may be derived for nonlocal condensates. As a specific example, the leading- 



D , order equations for the nonlocal condensates appearing in the quark propagator are derived 



and explicit solutions are obtained. Some applications of our analytical results are briefly 
discussed. 



PACS Indices: 12.38.Lg; 12.38.Cy; 12.38.Aw. 



I. Introduction 

The problem of strong interaction physics has been around for more than half a century, 
but the very nature of the problem varies with the so-called "underlying theory" which 
nowadays is taken universally to be quantum chromodynamics or QCD. Although the 
asymptotically free nature of QCD allows us to test the candidate theory at high energies, 
the nonperturbative feature dominates for hadrons or nuclei at low energies. As of today, 
it remains almost impossible to solve problems related to hadrons or nuclei, except perhaps 
through lattice simulations but it is unlikely that, through the present-day computation 
power which is still not quite adequate, simulations would give us most of what we wish to 
know. However, the need of solving QCD directly has gained tremendous importance over 
the last decade and, to this end, any fruitful attempt in confronting directly with QCD 
should be given due attention. 

The ground state, or the vacuum, of QCD is known to be nontrivial, in the sense 
that there are non-zero condensates, including gluon condensates, quark condensates, and 
perhaps infinitely many higher-order condensates. In such a theory, propagators, i.e. causal 
Green's functions, such as the quark propagator 

tSt^ix) =< I T{qt{x)q'^{0)) \ >, (1) 

carry all the difficulties inherent in the theory. Higher-order condensates, such as a four- 
quark condensate, 

< I T{i(;{zh^i;{z)qt{x)q'^{0)) \ >, 

with t/j^z) also labeling a quark field, represent an infinite series of unknowns unless some 
useful ways for reduction can be obtained. As the vacuum, | >, is highly nontrivial, 
there is little reason to expect that Wick's theorem (of factorization), as obtained for free 
quantum field theories, is still of validity. Thus, we must look for alternative methods in 
order to obtain useful results. 



It is known that the equations for Green functions up to a certain order usually involve 
Green functions of even higher order, thereby making such hierarchy of equations often 
useless in practice. In this paper, however, we wish to show that, provided that we may 
use the large Nc approximation to treat condensates of much higher order, there is in fact a 
natural way of setting up closed sets of differential equations which govern the inter-related 
Green functions to a given order. We consider this as an important accomplishment, both 
because we can always go over to the next level of sophistication in order to improve 
the approximation and because the large Nc expansion has been shown to yield desirable 
results for describing hadron physics. 

II. Leading-order Equations for Nonlocal Condensates 

In light of the nontrivial QCD vacuum, we begin by considering the feasibility of 
working directly with the various matrix elements such as the quark propagator of Eq. 
(1). Useful relations may be derived if we regard the equations for interacting fields [1], 

{t7^{d, + tg^A';,)-m}i; = 0; (2) 

d'-G;, - 2gr'^G%A: + <7V y7mV' = 0, (3) 

as the equations of motion for quantized interacting fields, subject to the standard rule 
for quantization that the equal-time (anti-)commutators among these quantized interating 
fields are identical to those among non-interacting quantized fields. As our basic example, 
we allow the operator {i'j^d^ — m} to act on the matrix element defined by Eq. (1) and 
obtain 



{t^^d, - mUtSf^ix) =td\x)d^%,+ < I T{{g—A;rQ{xmq',{0)) \ > 



2 

(4) 



We should always keep in mind that the QCD vacuum | > is a nontrivial ground state 
which is in general not annihilated by operating on it the annihilation operators. 



Eq. (4) can be solved by splitting the propagator into a singular, perturbative part 
and a nonperturbative part: 



iS^^Hx)=zSl^^''\x)+zSt^Hx), (5) 

where 

.5^f"^x) = |^e---z5if^p), (6a) 



with d = 7^a„ for a four- vector a„. The nonperturbative part then satisfies the equation: 



{i^f^d, - mUiSf^ix) =< I T{{g—A]:,^^q{xmq',m \ > . (7) 

This would be pretty much the end of the story unless we could find some way to proceed. 
We should note that Eq. (4) may also be derived by making use of, e.g., the path- integral 
formulation, and the issue of how to define renormalized composite operators, i.e. products 
of field operators, is by no means trivial (and fortunately we need not worry about such 
problem for the sake of this paper). 

As a useful benchmark, we note that, with the fixed-point gauge, 

a-^{x) = -^g;^x'' + --; (8) 

we may solve the nonperturbative part iS^Hx) as a power series in x^, 

+ ^<m-Gq> 5-'x^5,, + ■ ■ -, (9) 

The first term is the integration constant which defines the so-called "quark condensate" , 
while the mixed quark-gluon condensate appearing in the third term arises because of Eqs. 
(7) and (8). It is obvious that the series (9) is a short-distance expansion, which converges 
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for small enough x^. We note that is just the standard quark propagator cited in most 
papers in QCD sum rules [2]. 

The approach which we suggest here [3] is based upon two key elements, namely, the 
set of interacting field equations plus the rule of canonical quantization (for interacting 
fields). The equations which we obtain, such as Eq. (5), are much the same as the set 
of Schwinger-Dyson equations (for the matrix elements). An important aspect in our 
derivation is that the nontriviality of the vacuum | > is observed at every step — a 
central issue in relation to QCD. 

As another important exercise, we may split the gluon propagator into the singular, 
perturbative part and the nonperturbative part. The nonperturbative part is given by [3] 

g'<o\:G;ux)G::^{oy.\o> 

i:nm 

= -^ < g'^G'^ > {9iic.9up - Qiipgua) 

- -r^ < g^G^ > {x'^{gfj,a9iyp - 9fip9ua) - g^iaXiyXp + Qfj^pX^Xa - QufiX^Xo, + QiyaX^^Xp} 

+ 0{x^), (10) 
with 

< g^G^ >=<0 |: /G';:,(0)G'"'^"(0): | >, (11a) 

< g^G^ >=<0 |: (7'/"VG'«,(0)G'''«^(0)G'^,(0): | > . (116) 

Again, the first term in Eq. (10) is the integration constant for the differential equation 
satisfied by the gluon propagator. Note that inclusion of the gluon condensate < g^G^ > 
in Eq. (10) is standard but the triple gluon condensate < g^G^ > is a new entry required 
by the interacting field equation (3). Our approach indicates when condensates of entirely 
new types should be introduced as we try to perform operator-product expansions to higher 
dimensions. 

In what follows, we wish to focus on the quark propagator, as specified by Eqs. (1) 
and (4)-(7). We write 

iSt^ix) = S'^'id^.fix^) + ix,jg{x^)}, (12) 
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with x'^ = Xq — x^. f{x'^) and g{x'^) are what we refer to as "nonlocal condensates" in 
connection with the quark propagator. We note that 

<:qix)qi0):>=-12fix''), (13) 

which is the nonlocal quark condensate in the standard sense. To proceed further, we work 
only with the leading term in the fixed-point gauge and introduce 

<: {gG^M^mq^iO) ■■> 

^^'^^{{ifj.Xu - 'juX^)A{x'^) + ia^^B{x'^) + (7^x^ - 'y„Xf^)xC{x'^) + ia^„xD{x'^)}, (14) 

with G^i, = (A'^/2)G'^j^, an antisymmetric operator. The invariant functions A{x'^), B{x'^), 
C(x^), and D{x'^) are additional nonlocal condensates which we deal with explicitly in this 
paper. 

Under the assumption that we keep only the leading term in the fixed-point gauge (a 
simplifying assumption which can be removed whenever necessary), we have 



{i7-c>« - mU <: {gG^Mx)rkq'^{0) :> 
-x^^ <: ig^G^Mc,p-f'^q[x)l'^i,^, 



- - ^x^ <■■ {9'G^,G^0rqix)}1q^{O) :> (15a) 

1 1 4 



2 96 3 <^'^'><^{(7a*^--7.^m)?(^)}"^'(0):> (156) 

< g^G^ > d^\-f^x, - -f,x^){f{x^) + txg{x^)}. (15c) 



144 

Here the second line, (15a), follows from the field equation and the third line, (15b), is 
based on the large Nc approximation that the contribution in which G^i, and Gap do not 
couple to color-singlet is suppressed by a factor of l/N^- Thus, the factorization in the 
present case is justified to order l/N^, rather than just order l/Nc- We also note that we 
have used Eq. (10) to leading order in obtaining Eq. (15b), but this approximation can 
be relaxed if necessary. 

Now, we may use Eqs. (7) and (15c) and obtain a closed set of equations: 

2f'{x^)-mg{x^) = --i{B-x'^C), (16a) 
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2x^g\x^) + 4g{x^) + mf{x^) = ^x'^{A - D), (166) 

MB' - 2iC - 2ix^C' -mA = —^ < g^G'^ > /(x^), (16c) 

2iA + 2ix'^D' -niB ^0, {16d) 

- 2iA' + AiD' -mC = ---^< g^G^ > g{x^), (16e) 

2iB' + 2iC - mD = 0, (16/) 

where the derivatives are with respect to the variable x'^. 
Treating m as an expansion parameter, 

oo 

F{x^) = Y,m'Fk{x^), (17) 

fc=0 

we may solve the coupled equations, (16a)-(16f), order by order in m. To leading order in 
m, we obtain 

^'fo" + 3/o' - ^o^'/o - 2^o/o = 0, (18) 

{x'fg'^' + bix'fg'^ + {2x' - ^o'(^')'K " {2 + 2eo{x'f}go = 0, (19) 

with ^0 =< 9^G^ > /384. The equations for Aq, Bq, Cq, and Dq can easily be solved once 
we obtain /o and go- 

Eq. (19) can be simplified considerably by introducing 

go{x^) = (x2)-2^o(^'), (20a) 

which leads to the equation: 

x^~g'^' - ~g'i - iW~9o = 0. (206) 

Eqs. (17) and (20) can be solved by iteration, leading to the result: 



2~4^^'^^' ^ 2^ . 6 



/o(t) =ao{l + -r^i^ot)^ + -^^^^^{^,t)^ + ...} 

+ait{l + ^{Cot? + w-^iCot)' + ■■■}, (21) 



~g',{t) =c,e{l + ^(^ot)^ + jT^i^otf + •■■}, (22) 



with t = x^ and 

ao = -Y^<gg>, ai = Y^<q9c(T-Gq>, C2 = ^^-^4 — . (23) 

Eq. (23) is obtained by comparing to the weU-known series expansion for the quark 
propagator (see, e.g., [4]). Note that there are two integration constants, ao and ai, for 
f{x^) but there is only one permissible constant for g{x^) (and C2 is in fact a four-quark 
condensate taken in the large Nc limit). 

For a number of applications, it is useful to obtain analytic expressions for fo{t) and 
goit). This turns out to be possible by way of Laplace transforms. 

/■OO /"OO 

Us) = / rfse-^Vo(t), %{s) = / dse-'%(t). (24) 

Jo Jo 

We obtain 

fois) = -^ - ^^^sec-'^ + ^3^, (25) 

with 7o = 2ai/^o- Looking up the table for Laplace transforms, we find 
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9o = ^-Cot-Ii{Cot): (27) 

with /i (2) the modified Bessel function of the first kind, to order one. It is straightforward 
to show that Eq. (27) yields the series expansion in Eq. (22). Also, the function /i(^o^) 
enters the second series in fo{t) as in Eq. (21). 

To close our presentation of the explicit solution to leading order in m, we note that 

Eqs. (16) yields 

m) = -z^itBoit))' , (28a) 

Co(t) = -S^(t), (286) 

tg'oit) + 2go{t) = -^t{tDo{t)y, (28c) 

Aoit) = -tD'oit). (28rf) 



Thus, the functions Aq, Bq, Cq, and Dq can be solved exphcitly once /o and go are known. 
On the other hand, we may go beyond the leading order in m and obtain, as example. 



tf[" + 3/r - ^otn - 2^0/1 

= ltg'^ + ^90 - litA'o + 2Ao + t{tDo)" + 2(tI)o)'}; (29a) 

t^g'l' + btg'l + (2t - eGt^)9'i - (2 + 2^2^')^i 
= - 2 (^'/o' + tf'o - /o) + ^(-^^0 + 2^Co + ttC'^). (296) 

These equations can also be solved by Laplace transforms. In other words, the nonlocal 
condensate functions f{x^) and g^x^) can be analytically solved order by order in m. 



III. Discussion and Summary 

Thus far, we have described how to obtain a closed set of coupled equations for the 
nonlocal condensates which are relevant in the description of the quark propagator. We 
have also shown how these equations can be solved explicitly. Furthermore, some of the 
assumptions underlying our equations can be relaxed and more elaborate equations may 
then be obtained. Of course, some of our results are gauge dependent as the quark propa- 
gator (1) has been analyzed in a specific gauge (8). Nevertheless, our primary motivation 
for studying the quark propagator stems from our interest in the method of QCD sum 
rules [1], which may be regarded as the various approaches in which one tries to exploit 
the roles played by the quark and gluon condensates for problems involving hadrons. In 
this regard, results based on our leading-order equations are often adequate and the final 
answers are in general free from the potential gauge dependent problem. 

There are several approaches of QCD sum rules towards hadron physics. As the first 
approach, we may consider the Belyaev-Ioffe nucleon mass sum rules [5,4], where the short- 
distance expansion for the quark propagator is needed up to a certain (high) dimension. 
In this context, our analytical results on nonlocal condensates may not be very useful if 
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the resultant series converges rapidly, and in general our analytical expressions may be 
used to perform further analytical analysis of the problems as a way to improve the results 
obtained via short-distance expansions (in x^). The second approach is to consider the 
response of the QCD vacuum to some external fields, such as the method of QCD sum 
rules in the presence of an external axial field Z^{x) [6]. In this context, certain induced 
condensates are introduced (previously as new parameters) but the method offers a simple 
extension of the first approach in calculating magnetic moments, coupling constants, and 
other quantities by avoiding a need to treat explicitly the three-point Green's functions 
- a need which would still involve some conceptual difficulties. What is of great interest 
is that our analytical expressions for nonlocal condensates help to determine the induced 
condensates previously treated as new parameters, thereby making the external-field QCD 
sum rule method more powerful than what it used to be. To illustrate the point, we 
consider the external axial field Z^ with the interaction, 

SC{x) = gZ^{x)q{x)^^^^q{x). (30) 

For a constant Z^ field, there are two major induced condensates [6]: 

< I ^(0)7m75?(0) I >z« and < | q(0)gcG^.7"q(0) | >z« . 

We now have 

< I ^(0)7^75^(0) I >z« 
=1 fd^xgZ'^ix) < I T(Q(x)7,75Q(x)Q(0)7^75g(0)) | > 

=i J d^xgZ"{x){Tr[iS^'^\-xh^^,iS^''\xh^^,] 
+ Tr [z^(-x)7a 75^5''^°^ (x)7^ 75] 
+ Tr[z5'^°)(-x)7a75Z^(x)7^75] 

4-Tr[i^(-x)7a75«^(x)7^75]}. (31) 
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The first term is the one-loop result which can be regularized (e.g., in d dimensions) and, 
as expected for a perturbative contribution, its finite part is small compared to the second 
and third terms. The last term, which can be treated numerically, involves products of two 
condensates and it is higher than the second or third term by at least three dimensions 
(and is likely of less numerical significance). 
Analogously, we have 

< I m9cG^.Yq{0) I >zo^ 

=i fd^xgZ'^ix) < I T{q{x)-f^-f5qix)qiO)gcG^,YqiO)) I > 
^igZ'^g^^e^^x^ J d^xTr{^^^,iS'-ix)Y <: {G^^^ q{0)r q\x) :>}. (32) 

Thus, our analytical results on Aq, Bq, Cq, and Dq [Cf. Eqs. (14) and (28)] may be used 
to evaluate Eq. (32). 

Without going into the details (which shall be presented elsewhere together with 
detailed discussions), we mention that, upon Wick's rotation on the time integration 
(J_ ^ /_ ), our analytical solutions can be employed explicitly. The final results 
for the above two induced condensates are recorded immediately below: 



< I ^(0)7m75?(0) I > = gxZ^ < qq >, (33a) 

< I m9cG^r.rq{0) \0>= gnZ^ < qq >, (336) 

(^\ w? ,, Tim? 

;^(i) < gg > = In ^ + 7s , 33c 

(2) _ TT m _ 1 mml gl <qq> 

"" "2^0 16 eo ^216 eo ' ^ ^ 

^ = ^(.r + 21n2 + l)--L^. (33e) 

Numerically, we find xa ~ 0.15 Gey^ and Ka ~ 6 x lQ~^GeV'^. (Here we have used 
the updated value on the gluon condensate [7], the current quark mass, and other input 
parameters as adopted previously [4].) Such values for the induced condensates, albeit 
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somewhat smaller than the commonly adopted ones, are not unexpected since the leading 
contributions are linear in the "current" quark mass m. 

There are other problems for which our analytical results may be very useful. For 
example, the amplitude given by 

T^Aq^p-q)=i fd^xe-^'^-^ < n+{p) \ T{J^{x)MO)) \ n+{p) > (34a) 

= {-9fiiy + -^)Ti{q^,p-q) 

1 T) ' CI T) ' (1 

+ — (Pm T^i^)^P^ ^qiy)T2{q'^,p-q), (346) 

characterizes the forward Compton scattering off tt"*" (a Goldstone boson) and also the 
parton distributions of tt"*". Applying the soft-pion theorem (together with current algebra), 
we find, in the limit of p^ ^ 0, 

T^^Aq^O) 

jd'xe-^'^-^ < I T{{Al{x) - tAl{x)}{Al{0) + tA^O)} - 2V^{x)V^\0)) \ > . 



f2 

J IT 

i 

'■J2 



d^xe-''^-^Tr{iS'/{-xh^T,tS:\xh,r 



- -z5^«(-x)7^z5r(^)7. - -tST{-x)-i,iSf{x)-f,}. (35) 



The structure functions Wi{q^,p.q) (in the description of deep inelastic scattering off the 



TX 



+ target) is the imaginary part of Ti{q^ ,p • q) divided by the factor tx. Eq. (35) suggest 



that our analytical expressions for the nonlocal condensates may be useful for analyzing 
properties of Goldstone pions. We find, as a soft-pion limit, 

Wi{q^,p-q) = -ImTi(q2,p-q) 

TT 

— ' of2t ' ^^ P/x ^ and q^ -^ 0. (36) 

This limit is derived making use of our analytical expressions on the nonlocal condensates, 
again with Wick's rotation on the time integration. 
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In summary, we have in this paper suggested a specific way of obtaining a closed set 
of coupled differential equations for nonlocal condensates. Specifically, the leading-order 
equations for the nonlocal condensates in relation to the quark propagator are obtained 
and explicit analytical solutions are obtained. We believe that such results could be very 
useful for a large number of problems in hadron physics. 
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